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Abstract 

We give a simple presentation of the combinatorics of renormalization in pertur- 
bative quantum field theory in terms of triangular matrices. The prescription, that 
may be of calculational value, is derived from first principles, to wit, the "Birkhoff 
decomposition" in the Hopf-algebraic description of renormalization by Connes and 
Kreimer. 
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1 Introduction 



A Hopf algebra structure underlying the combinatorics of perturbative renormalization was 
recognized by Kreimer [1]. Some have worried about the practical usefulness of his insight 
for organizing everyday computations in quantum field theories. A partial answer to this 
legitimate question is given in [2] , for instance. This paper presents another partial answer, 
in a different vein. We show that the Feynman rules collectively possess a triangular matrix 
representation, such that the renormalization map becomes a matrix operation. This infinite 
matrix can be truncated almost ad libitum. 

The procedure is general, largely independent of the renormalization scheme (although 
we illustrate everything with the MS-scheme in dimensional regularization), and essentially 
independent of the particular field theory model one works with. The latter enters just in 
filling up the matrix entries, generally entailing further simplification. Only knowledge of 
linear algebra and quantum fields are required; no Hopf algebra is (openly) used. More detail 
is found in [3]. 

First we recall the algebraic behaviour of the subtraction map K: the whole paper 
turns around its Rota-Baxter property. In Section 3 i^' is lifted to the matrix level, and 
the computational recipe for the matrix counterterm and matrix renormalization maps is 
found. In the following section we verify that this reproduces the diagrammatic Bogoliubov 
operation; we take examples from the 0| model, and compare with the tables in [4]. In 
Section 5 we rework the matrix representation using the map that picks out the finite 
parts. The next two sections contain mathematical summaries. In Section 7 we sketch our 
derivation of the matrix representation. Finally we examine the outlook. 



2 The subtraction map as a Rota— Baxter operator 

Consider Laurent series 

Sie) = hao + OieH . (1) 

With the ordinary multiplication, they form a commutative algebra V with unit. Consider 
further the operation K which picks out the pure pole part 

mie) = — + -;;^ + --- + — 
and the operation :— id — K keeping the finite part, 

K_^.[S]{e) = ao + aie-\ . 

The projector condition = K ensures that the intersection between K{y) and Kj^{y) 
is zero. The product of two elements of KiV) remains in KiV) — and likewise for Kj^iV). 
The key property 

K[S^] K[S2] = K [K[S^]S2 + S^K[S2] - ^1^2] , (2) 
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is easy to check — see Section 6. It makes K a Rota-Baxter operator [5]; also is a Rota- 
Baxter operator. All this applies in particular to series corresponding to dimensionally 
regularized integrals in the MS-scheme subtraction. Our arguments are purely combinato- 
rial, so we need not worry about the precise form of the coefficients. We adopt for the 
subtraction operator K the notation of [6], followed in [4]; sometimes we write for clarity. 



3 Setting up the recipe 



In this section we suppose given the pair {V, K) of a commutative algebra with unit and a 
Rota-Baxter projector, first as an abstract framework; we always have in mind the algebra 
of unrenormalized Feynman amplitudes of the form (0). Consider upper triangular matrices 
of finite size with entries in V , of two types: nilpotent, i.e., with O's on the diagonal, and 
unipotent, i.e., with I's on the diagonal: 



/O * 




Vo ... 



Zij 



0/ 



/I * 

1 



Vo ... 



l) 



We define Rota-Baxter operations /C_ = /C and /C+, on algebras M"p(V^) of upper triangular 
matrices with scalar diagonals and with entries in V , by extending the maps K and Kj^ 
comp onent wise , 

Verification for /C, /C+ of the analogue of is immediate; but the algebras M"P(y) are no 
longer commutative. We seek to factorize an arbitrary unipotent element if in the form 



ip = ip+Lp_ 



(3) 



where the factors ip- G 1 + /C[M^p(V^)], ip^ G /C+[M^p(V)] are also unipotent; note that they 
are unique. This can be called a matrix Birkhoff decomposition. 

If /C[log(y9] and /C+[log(y9] happened to commute, it would be enough to choose (p+ = 
g'c+[iog(/j] g^j^^ _ g-/c[iog</3]^ jj-^ general, that is not so; but we are able to compensate for the 

lack of commutativity between the images of /C_ and /C+. For that, consider the equations 

= 1 — /C_ [{if — l)ip-] and ip^'^ = 1 — /C+ — 1 

respectively solved by 



= 1 - /C_[v9 - 1] + /C_[(<^ - 1)/C_[v9 - 1]] - 
= 1 - /C+[v9 - 1] + /C+[/C+[^ 1)] - 



(4) 
(5) 



3 



Both series terminate. Atkinson's theorem [7] asserts that these matrices ip_^} verify (jH)). 
The proof runs as follows: 

^;V- = i-ic+ [^^\^ - 1)] - /c_ [iv - i)^_] + /c+ [^-\^ - i)]/c_ [iv - i)^_] 
= 1 - v+^v - 

after some work with the Rota-Baxter property. 

The matrix (f^ is what we are really after. It can obviously be obtained as ff^, or by 
inverting (jSJ, using the geometric series formula 

(^+ = i-(v9;i-i) + (^;i-i)2-.... 

A better course is perhaps to observe that, by the same token as in we obtain 

(p+ = l- /C+[((^"^ - l)ip+]. 
Thus, the respective formulas for the components of and are 

j—i—l 

= -K4^.,] E H'-'K- K ^-bw. • • • K.i^iJ ■■■]], (6) 

k=l i<li<h<—<lk<j 
j—i—l 

= -K^ W] + E E (-)'''^+ [^ii ...K^ brj-] •••]]• 

k=l i<h<h< - <lk<j 

These similar formulas are our workhorses; with the appropriate definition of if, the ma- 
trix will be seen to contain all the information on counterterms in renormalization; 
and ip+ on the renormalized quantities. 

4 Making the recipe work 

Now we make explicit how the Feynman rules specify such operators ip. Recall that if Fj C Fj 
is a superficially divergent subgraph of Fj, the cograph Fj/Fj is obtained by shrinking Fj to 
a vertex within Tj. We only consider subgraphs that are generalized vertices [6]. Chosen 
an ra-point function, the spaces of vectors on which the matrices act are spanned by the 
corresponding (superficially divergent, connected, amputated) Feynman graphs. We may use 
the familiar bra-ket notation to denote the diagrams as vectors. A basis |Fi), IF2), IF3), . . . 
for such a space can be ordered in many ways, the only conditions being that |Fi) = |0) 
— the empty diagram — and that each cograph of any F^ occurs in the basis as some F^ with 
m < I. It is then convenient to order the basis by number of loops (or vertices, if we work 
on coordinate space); but the order within a given loop- number sector is immaterial. Once 
the external structure and the basis are fixed, we fill up the entries of a matrix by the rule: 
for i ^ j, 

ifij = (unrenormalized) amplitude of F' if Fj ~ Fj/F' , 
r' 
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otherwise ^pij = 0. This entails triangularity, since (pij = if i > j. We set (pu = 1 for all i. 
Note that T' need not belong to the basis list (it might be disconnected, for one thing). Let 
<f{T') be the unrenormalized amplitude of T'. We just said that the coefficient of iFj) in 
ipdTj)) is '^j-i 'f>{T') for Ti ~ Tj/r'. The notation is appropriate because (p is the abstract 
object represented by the matrix ip (see Section 7); it has the property that 

<p(TiTj) := p(Ti U Tj) = (p(Ti) ipiTj). 

This property is shared by the elements of and (/?+. As a bonus, it allows to simplify the 
notation later on, by simply omitting p. 

When we truncate the matrices we are not obliged to include all the diagrams belonging 
to the higher sector — and we can also choose, for whatever purpose, particular classes of 
diagrams, subject to the aforementioned two conditions. With that, the first row of p is 
given by 1, (^(r2), (^(Fs), . . . , the unrenormalized amplitudes of all the diagrams; and the 
analogously defined first row (1, (^„(r2), (^-(Fs), . . . ) of will yield all the counterterms of 
the theory! 

We take as a simple example the space of graphs relevant to the 4-point function for the 
0| model, truncated to the 12 tadpole-free diagrams (including the empty one) up to three 
loops. We adopt the order of [4] for the basis, as follows: 

Ti r2 r4 F5 Fg F7 Fg Fg Fio Fn F12 

O =00 y 000 <M5 xgx ^ O ><II> 

Notice that F2 = F7/F' where the sunset diagram F' = -Q- does not appear in the basis list. 
We find that p is equal to 

/l <f{T2) (fiTs) ifir^) ^(Fg) (^(Fg) ifiT-i) ip{Ts) (^(Fg) ^(Fio) (^(Fn) (^(ri2)\ 

1 20{T2) ^(Fa) 2^{r3) + ^HT2) ^{T^) + ^^T2) ^{T') ^(Ts) ^HT2) ^{T^) 2^(r4) 
1 3^(r2) (^(Fs) f>iT2) 

1 (^(Fs) 2^(r2) 2f>{T2) (^(Fa) 

1 000000 

1 

1 

1 

10 

10 
1 

V 1 / 

The abundance of zeros is welcome in the calculation: with this concrete form of p the 
series (0} stops after three iterations. Omitting p from the notation as advertised, one reads 
off from (jU)) the first row of containing the counterterms for the eleven nontrivial dia- 
grams, that we write as the column matrix f'^ji in the following display. These expressions 
for (^_(Fj) coincide with those listed in the tables in [4], where they are denoted —KR(Tj). 
But note that here the Bogoliubov preparation map R does not appear explicitly; our re- 
sult, with one exception, is not recursively presented, and, with the help of some symbolic 
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programming, can be obtained at one stroke. It is clear that the method will jointly handle 
large numbers of multiloop diagrams with ease. (Of course, we are not claiming that it is 
always quicker than the standard procedures.) 



/ 1 \ 

^-(o<0) 

^_(><gx) 



( 



1 

-ir[Y] + ir[>Oi^[>0]] 
-i^[>0]^-(Y) 



\ 



V-^[>^] + K[<y'K\<y\\ + 2K[)Oi^[ Y]] - 2[ir[>o]]^3V 



Above we wrote [-ftr[)0]]''^^^ := -ft'fO -^[>0 -^[O]]]- The graphs r3,r5,r6 are cutvertex, 
for which (p and the renormalization map are known to factorize. For them (jH)) prima facie 
gives a more complicated expression, that can be reduced to the expressions shown by some 
Rota-Baxter gymnastics. We give the example of (O<0 ) . From our matrix operations: 

^_(>o<0) =-K[<x^]+K[YK[<y]\+K[<yK[Y]] 

Since K[0^] = K[<yY] = -K[>0]K[Y] + K[K[>0]Y] + ^[^^ [ Y] 'O] , we get 

<^_(0^) =K[^]K[Y] +i^[0'i^[0]] +K[^K[Ky^]\ -3[ir[>0]]^'^ 
= K[o\K[Y] + K[<y'] + K[<y^] K[o] - 3[K[<y]]^'l 

To continue, we invoke the classical Bohnenblust-Spitzer identity [8]: 



n\ [K[A]] 



K[AK[A...K[A] 



n times 



Pe u„ pe P 



which is itself derivable from the Rota-Baxter identity; here n„ is the set of partitions p of 
the set {1, . . . , n}. For the present case. 



6[ir[0<]]^^^ = K[>cy] K[K>] K[>0] + 3K[Ky^] K[Ky] + 2K [o<^] , 



implying that 

-K[<y]K[KyK[Ky]] = {-k[k>]) {-k[Y] + K[K>K[<y]]) =^-{<y) ^-(Y). 

According to the theory underlying the matrix representation (Section 7), the factorization 
property of if- is automatic: we did not prove anything, just performed an internal check. 



5 The matrix representation in terms of K+ 

Renormalization theory is usually formulated in terms of subtractions [4,6,9]. For good 
reasons: for instance, in the MS-scheme the counterterms are local, independent of the mass 
and the renormalization scale; this helps to establish renormalization group equations. Also, 
the derivation of the whole procedure is simpler in terms of subtractions. However, the 
understanding of renormalization as an approximation process, rather than a cancellation of 
infinities, is thereby lost. Recently 't Hooft expressed the desideratum of a renormalization 
scheme exclusively containing dressed vertices [10]. We take here a step in this direction 
by rephrasing the renormalization of (regularized) Feynman graphs with sub divergences in 
terms of K^. The magic of the triangular matrix representation implies that the matrix cp^ 
must give, graph by graph, the completely renormahzed expressions {id— K)R(T) = K+R(T). 
Inspection of the equations in Sectional on the other hand, shows that the same calculation 
method for in terms of K_ yields (1, ip^(T2), (p+(Ts), ■ ■ ■) in terms of K+, provided one 
starts by inverting the matrix if. This is a small price to pay, and again a bit of symbolic 
programming goes a long way. We illustrate the procedure with the same model example: 



/ v^+(oo) \ 

MY) 

<P+ (oco) 
^+(00) 

^+(x2x) 
MU) 
^+ (<(]>) 



V 



/ K4Ky]K4<y] 

k4Y]+k+ [<y] \o] - \OK^ \o]] 

K+[K>]K+[Ky\K+[Ky\ 
[O] [ Y] + K40] ~ K40K40]]) 

^+[><g><] +^+[o<]i^+[0]-^+[0^+[o<]] 

_ [^ + 2ir+ [ Y ] [O] - [O] [O'] 

- 2ir+[>Oi^+[ Y]] - K4o^K4k>]] + 2[ir+[0]i^'^ 
_ [O] + 2i^+ [ Y] KAo] - KAO] KAO'] 

- 2;^+[o<i^+[ Y]] - k4k>'k4k>]] + 2[ir+[o]i^'^ 

K+ [><D>] + 2ir+ [ Y ] ^+ [o] - [o<] [o' 
-2K^[y:K4<y]]+K^[<y{K4<y]f] 



\ 



K 
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The simpler cases have been omitted. The reader is reminded that -ft'+fOx] actually means 
-^+[<^(0<)], and so on. 

6 More on Rota— Baxter operators 

A few extra comments on the Rota-Baxter property of and Kj^ are in order. In the work 
by Kreimer, the Rota-Baxter property appears for the first time in [11], under the name 
"multiplicativity constraint". The maps K^^ may be regarded as generalized integrals. 
Indeed, let us insert a parameter 6 (a Rota-Baxter weight) before the last term of (P)): 

K[Si] K[S2] = K[K[Si]S2 + SiK[S2] - eSiS2] . 

The case 9 = corresponds to a property of the integral I[f]{x) := f{t) dt, to wit, 

/[/l]/[/2]=/[/[/l]/2 + /l/[/2]], (7) 

which is just integration by parts. For g fixed, the solution of the equation / = 1 — I[gf] 
with / satisfying ^ is given by 

f=l-I[g]+I[gI[g]]-... = e-'^^\ 

which follows from (|7j) for /i = /2 = /, and illustrates our approach in this paper. 

Although (j21) is elementary, we wish to prove it here. This is warranted because the 
Rota-Baxter property is persistently ignored in field theory treatises; and this neglect is 
not without consequences. For instance, in section 5.3.3 of the standard text [9], we find a 
tortured argument to try to prove (^-(Fj U F^) = ip_(Ti)ip_(Tj), in which the intermediate 
formulas (5.3.15) and (5.3.16) are plain wrong. To see why (0) holds, notice that 

K[Si]S2 + SiK[S2] - S1S2 = K[Si]iK[S2] + i^+N) - iK[Si] + K4Si])K4S2] 

= K[S^]K[S2]-K4S^]K4S2], 

and applying K to this equality kills the term _ft'+[S'i] K^[S2], leaving -ft'[S'i] K[S2] unchanged. 

7 The rationale for the matrix representation 

The not so mathematically-minded might wish to skip this section. In the formalism Kreimer 
developed jointly with Connes [12, 13], Feynman diagrams are organized in a Hopf algebra 
Tijf of graphs; Feynman rules are understood as linear and multiplicative maps (p oi Tijr 
into an algebra V (commutative, with unit) of quantum amplitudes; and the disentangling 
of sub divergences is formulated as a factorization problem (Birkhoff decomposition). The 
original version had a strong geometrical flavour, but its supporting algebraic frame has 
emerged since [5]. 

The space Tijr is the algebra of polynomials with connected Feynman graphs as indeter- 
minates, multiplication being given by simple juxtaposition of graphs. Connes and Kreimer 
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introduced on Tij: a coproduct A : Tij^ — > Tijr ® Hj^, serving to encode the superficially 
divergent subgraphs, by setting A(r) := J^r ^' ® T/F', in the notation of Section El For the 
coproducts of a Hopf algebra H one writes Aa = J2 "^(i) ® '^(2); for a ^H. 

Let {V, K) be a commutative Rota-Baxter algebra, and consider Hom(H, V), the space of 
linear maps from Ti to V] this is an algebra with the convolution operation, given by f-kg[a) = 
^ /(a(i)) (7(a(2)), for a E H. In our case the multiplicative (that is, product-respecting) 
elements of Hom(7ijc-, y), with V the algebra of Feynman amplitudes, are of particular 
interest. Clearly they are determined by their action on the subspace JF of connected graphs. 
We construct a representation \1/ of Hom(?-^jF, V) by infinite triangular matrices with entries 
in V by taking the composition 

where my is just multiplication on V. The plot works because the external structure of the 
cographs F/F' is the same as that of F, so A actually sends J-' into TCj^ ® J-". Thus for any 
/ G Hom(?ijc-, V) a connected graph is sent by into a linear combination of connected 
graphs with coefficients in V, corresponding to the same n-point function. In fact, is an 
antirepresentation, since -k g] = 

With the operator K given by K[/](a) := K[f{a)], the space Hom(HjF, "K) becomes a 
(noncommutative) Rota-Baxter algebra; then \1/[K[/]] = /C[\&[/]], with /C the known matrix 
Rota-Baxter map. Let finally ip G Hom(7ijr, V) be the Feynman rule, which is multiplicative, 
and denote 

if := ^[(^]. 

This will be a unipotent matrix. We have at last reproduced the setting of this paper. 
The matrix decomposition Q is a consequence of Connes' and Kreimer's algebraic Birkhoff 
decomposition ip = <^I^ *<^+, where the two factors are multiplicative as well [12]. Proofs 
and details are found in [3]. 



8 Conclusion 



Inspired by the Connes-Kreimer Hopf algebra formalism, we have exhibited the combina- 
torics of renormalization as a collective process, mechanized by means of simple matrix 
calculus. Our approach neatly resolves the tension between the "additive" and the "multi- 
plicative" sides of renormalization: the recursive diagrammatic subtraction of sub divergences 
is the outcome of a multiplicative process (this is not quite a trivial remark: a direct check 
of the relation <^+(F) = K^R[r) if and only if <^-(F) = —K-R{r) involves somewhat messy 
calculations). As a consequence, the renormalization of the Lagrangian's parameters by 
counterterms takes place by composition of series; the latter has been known since 1855 to 
have a triangular matrix representation [14]. All this is more or less clear from the analysis 
in [6,13,15]; but probably deserves further elucidation. 

Also, we have rewritten the renormalization map in terms of the projection on the 
finite part. In dimensional regularization this prescription falls short of 't Hooft's desider- 
atum [10], as some of the terms in contain coefficients of the pole parts; this objection. 
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nevertheless, loses force in regularization-free schemes like BPHZ and Epstein-Glaser renor- 
malization, that also possess a Rota-Baxter property [16]. 
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